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Two flavor color superconductivity and compact stars 
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Baryonic matter at high density and low temperature is a color superconductor. In real world, this state ol matter may 
naturally appear inside compact stars. A construction of a hybrid compact star with two flavor color superconducting 
quark matter in its interior is presented. 



1 Introduction 

At large baryon density, quantum chromodynamics 
(QCD) becomes a weakly interacting theory of quarks 
and gluons [ ^ . Because of an attractive interaction 
between quarks, the ground state of such quark matter 
is a color superconductor [I2IH1B]- At asymptotic den- 
sities, quark matter was studied from first principles in 
Refs. [0311)1 El Unfortunately, these studies are not 
very reliable quantitatively at realistic densities that 
exist in nature (i.e., at densities less than about lOpo 
where po « 0.15 fin -3 is the normal nuclear density). 

In general, QCD at high baryon density has a very 
rich phase structure. There are many possible color 
superconducting phases of quark matter made of one, 
two and three lightest quark flavors. Each of them 
is characterised by a unique symmetry breaking pat- 
tern and by a specific number of bosonic as well as 
fcrmionic gapless modes. 

In the rest of this paper, we are going to concentrate 
almost exclusively on matter with two quark flavors. 
The corresponding ground state of matter is the so- 
called two-flavor color superconductor (2SC). In this 
phase, the color gauge group SU(3) C is broken by the 
Andcrson-Higgs mechanism down to SU(2) C subgroup. 
With the conventional choice of the condensate point- 
ing in the "blue" direction, one finds that the conden- 
sate consists of red up (u r ) and green down (d g ), as 
well as green up (u g ) and red down (d r ) quark Cooper 
pairs. The other two quarks (ub and db) do not par- 
ticipate in pairing. This is the conventional picture of 
the 2SC phase [CUHE]. 

In passing, we mention that quark matter at high den- 
sity may also contain strange quarks. This would 
be the case when the constituent medium modified 
mass of the strange quark is smaller than the value 
of the strange chemical potential. The current limited 
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knowledge of QCD properties at finite density does not 
allow us to resolve the issue regarding the strangeness 
content in baryonic matter at densities existing inside 
compact stars unambiguously. As a benchmark test, 
here we study only non-strange quark matter. 

We start our discussion by pointing that matter in the 
bulk of a compact star should be neutral with respect 
to electrical as well as color charges. Also, such mat- 
ter should remain in /3-equilibrium. Satisfying these 
requirements impose nontrivial relations between the 
chemical potentials of different quarks [E0E3- ^ n turn, 
such relations influence the pairing dynamics between 
quarks, for instance, by suppressing conventional 2SC 
phase and favoring the so-called gapless 2SC (g2SC) 
phase [ITT]. 

The condition of charge neutrality should not neces- 
sarily be satisfied locally. It is acceptable, for example, 
if matter stays in a mixed phase that is neutral glob- 
ally, or on average [E1EJ- Below we use this idea to 
obtain a model equation of state of non-strange hybrid 
matter that is made in part of the 2SC phase [I14|. 

2 Quark model 

We use the simplest SU(2) Nambu-Jona-Lasinio model 
of Ref. [El to describe quark matter, 



C = q(i^d^ - m )q + G S [{qqf + (^75^) 2 ] 
+ G D [(ifee b l5 q)(iqee b l5 q c )] , 



(1) 



where q c = Cq T is the charge-conjugate spinor and 
C — i7 2 7° is the charge conjugation matrix. The 
quark field is a 4-component Dirac spinor that car- 
ries flavor (i — 1,2) and color (a — 1,2,3) indices. 
t = (t 1 ,-!- 2 ,-!- 3 ) are Pauli matrices in the flavor space, 



and (e) 



k — 



b\af) — e af3b 



are antisymmetric 



tensors in flavor and color, respectively. We intro- 
duce two independent coupling constants in the quark- 
antiquark and diquark channels, Gs and Go- Also we 
restrict ourselves only to the chiral limit (mo = 0). 
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The values of the parameters in the NJL model are 
the same as in Ref. [ 16 : Gs — 5.016 GeV~ 2 and 
the cut-off A = 653 MeV. The strength of the diquark 
coupling Gd is taken to be proportional to the quark- 
antiquark coupling constant, i.e., Gd — i]Gs with r\ = 
0.75. The choice r\ — 0.75 corresponds to the regime 
of intermediate strength of the diquark coupling [ lllj . 

In /3-equilibrium, the diagonal matrix of quark chem- 
ical potentials is given in terms of baryon (fj,B = 3/z), 
electrical and color chemical potentials as follows: 

2 

A%>/3 = — ^eQij)5 a p H — 7^^sSij(Ts) a p, (2) 



The solution that satisfies the gap equation and both 
neutrality conditions was studied in detail in Ref. [ lllj . 
It was found that the corresponding phase of matter is 
the g2SC phase. This phase has the same symmetry of 
the ground state as the conventional 2SC phase. In the 
low-energy spectrum, however, it has two additional 
gapless fermionic quasiparticles. 

Our analysis shows that the neutral g2SC phase is 
more favorable than the neutral normal phase of quark 
matter [lllj. We find, in particular, that the pressure 
difference of these two phases is of order 1 MeV/fm 3 . 

3 Mixed phases 



where Q and Tg are the generators of U(l) em of elec- 
tromagnetism and U(l)s subgroup of SU(3) C group. 

In the mean field approximation, the effective po- 
tential for zero temperature quark matter in 0- 
cquilibrium with electrons takes the form [ ITT] 



n = n 



I'c 



12tt 2 4Gc 



4G 



D 



E 



d 3 p 



£«|,(3) 



where is a constant added to make the pressure of 
the vacuum vanishing. The sum in the last term runs 
over all (6 quark and 6 antiquark) quasiparticles. The 
dispersion relations and the degeneracy factors of the 
quasiparticles read 



E± = E(p)±n ub , 
E± = E(p)±n db , 



[xl] 

[xl] 
[x2] 



E± ± = v/[E(p)±/i] 2 + A2±^. 
Here we used the following shorthand notation: 



(4) 
(5) 
(6) 



E(p) 
ft 



\/p 2 + TO 2 , 



/U8 



f-^dg f^ur f-^dr f-^ug 



fie 
2 



(7) 
(8) 

(9) 



The value of the thermodynamic potential that deter- 
mines the pressure, fiphys = — P, is obtained from 
in Eq. J3J after substituting the order parameter A 
that solves the gap equation, i.e., dil/dA = 0. 

Now, if we restrict ourselves to the case of homoge- 
neous quark matter, the following two conditions of 
charge neutrality should be imposed: 



n 8 



dfi 8 



»Q 



dfi e 



(10) 



In addition to homogeneous phases, one could also 
study various mixed phases of quark matter. The most 
promising components for constructing mixed phases 
are (i) normal phase, (ii) g2SC phase, and (iii) ordi- 
nary 2SC phase. Note that the first two of them allow 
locally neutral phases by themselves, while the last 
one tends to be positively charged. In principle, other 
phases could also be taken into consideration. 

Inside mixed phases, the charge neutrality is satisfied 
"on average" rather than locally. This means that dif- 
ferent components of a mixed phase may have non-zero 
densities of conserved charges, but the total charge of 
all components vanishes. In absence of local neutral- 
ity, the pressure of each phase could be considered as 
a function of baryon chemical potential, as well as a 
function of chemical potentials related to other con- 
served charges (e.g., fx e and /z 8 in the model at hand). 
This information is used as an input to obtain the 
most favorable mixed phase construction (see below). 

While we intend to relax the local neutrality condition 
with respect to the electrical charge, in this paper we 
always enforce the condition of local color neutrality. 
Here we cannot fully justify this constrain. However, 
one could speculate that separation of color charges is 
less likely in mixed phases. 

Without going into much detail, let us give a brief 
introduction into the general method of constructing 
mixed phases by imposing the Gibbs conditions of 
equilibrium [ I12| which are the conditions of mechan- 
ical and chemical equilibrium between different com- 
ponents. In regard to a mixed phase of the normal and 
2SC components of quark matter, for example, these 
conditions read 



Mi Me 



(NQ) 



(2SC) 
1 i 

{ISC) 



(11) 

(12) 
(13) 



These conditions are easy to visualize by plotting the 
pressure as a function of chemical potentials (fj, and 
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Figure 1. Pressure as a function of /i = /is/3 and fi e for 
the hadronic phase (at the bottom), for the 2SC phase (on 
the right in front) and the normal phase of quark matter 
(on the left). The dark solid line follows the neutrality line 
in hadronic matter, and two mixed phases: (i) a mixture of 
the hadronic and normal quark phases; and (ii) a mixture 
of the normal and 2SC quark phases. 



/x e ) for each component, see Fig. ^ As is clear, the 
Gibbs conditions are satisfied automatically along the 
intersection lines of different pressure surfaces. 

In a mixed phase, different components occupy differ- 
ent volumes of space. To describe this quantitatively, 
we introduce the volume fraction of the normal phase 



of quark matter: X2SC 



Vnq/V (notation Xb 



"volume fraction of phase A in a mixture with phase 
B"). Then, the volume fraction of the 2SC phase is 
given by X 2 nq = (1 — X-zsc)- F rom t ne definition, it is 
clear that < X2.sc — 1- 

The average electric charge density of the mixed phase 
is determined by the charge densities of its components 
taken in the proportion of the corresponding volume 
fractions. Thus, 



4 MP) = X^cni NQ) 



+ (1 - X» 



(14) 



If the charge densities of the two components have 
opposite signs, the global charge neutrality condition, 
n ( MP ) _ can j-^ i m p 0sec | p Otherwise, a neutral 
mixed phase could not exist. In the case of quark 
matter, the charge density of the normal quark phase 
is negative, while the charge density of the 2SC phase 
is positive along the line of the Gibbs construction 
(dark solid line in Fig.^). So, a neutral mixed phase 



exists, and the volume fractions of its components are 



X2SC — 



V 25C 

Xnq 



(2SC) (7VQ)' 
He 1 be 



NQ _ n e 

2SC ~ ANQ) „(2SC) ■ 
lie i^e 



(15) 
(16) 



By making use of these expressions, we can calculate 
the energy density of the corresponding mixed phase, 



,(MP) _ NQ (NQ) 



NQ \J2SC) 



(17) 



This is essentially all that we need in order to con- 
struct the equation of state of the mixed phase. 

So far, we neglected the effects of the Coulomb forces 
and the surface tension between different components 
of the mixed phase. In a real system, however, these 
are important. In particular, the balance between the 
Coulomb forces and the surface tension determines 
the size and geometry of different components of the 
mixed phase. In our case, nearly equal volume frac- 
tions of the two quark phases are likely to form al- 
ternating layers (slabs) [ E| • The thickness a of the 
layers scales as <7 1 / 3 (ni 2 ' SC ' ) — ni N ^)~ 2 / 3 where a is 
the surface tension. In the model at hand a « 10 fm. 

The energy cost per unit volume to produce the layers 
scales as a 2 / 3 (n? SC) - n{ NQ] ) 2 ' 3 [|T7j. Therefore, the 
mixed phase is favorable only if the surface tension 
is not too large. Our estimates show that a max < 
20 McV/fm 2 . For slightly larger values, 20 < a < 
50 MeV/fm 2 , the mixed phase is also possible, but 
its first appearance occurs at higher densities, 3po < 
Pb 5po- Note that the value of the maximum surface 
tension is of the same order as in Refs. [ 1171 118| . In 
this study we assume that actual value of the surface 
tension is not very large. 

The validity of the quark model is limited when the 
baryon density decreases. In fact, at sufficiently low 
densities, quarks are confined inside hadrons, and it is 
more natural to use a hadronic model. Here we use 
the chiral SU(3) L xSU(3) i? , model of Ref. [HHj. 

It is expected that the hadronic and quark phases are 
separated by a first order phase transition. Then, the 
hadronic and quark phases can co-exist in a mixed 
phase [EI that is constructed by satisfying the Gibbs 
conditions similar to those discussed earlier. Thus, we 
plot the hadronic surface of the pressure along with 
the quark surfaces in Fig. 2] Again, the intersection 
lines of different surfaces indicate all potentially viable 
mixed phases. Although the Gibbs conditions are sat- 
isfied along all such lines, not all of them can produce 
neutral phases (e.g., there are no neutral constructions 
along the lightly shaded solid lines in Fig. [TJ. 
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The dark solid line in Fig. Ogives the full construction 
that consists of three pieces. The lowest part of the 
curve (up to the point denoted by □) corresponds to 
the neutral hadronic phase. Within this region matter 
is composed mostly of neutrons and a small fraction 
of protons and electrons. 

The mixed phase of hadronic and normal quark mat- 
ter starts at the baryon density ps ~ 1.49,00 (D-point 
in Fig. At this point the first bubbles of decon- 
fined quark matter appear in the system. At the be- 
ginning of this hadron-quark mixed phase, the decon- 
fincd bubbles are small but highly negatively charged, 
whereas the hadronic phase, in which the bubbles are 
embedded, is slightly positively charged. The charge 
neutrality condition reads 

ee xl Q n^ + (1 - x N H Q )n^ = 0. (18) 

where ni and He ^ are the charge densities of 
hadronic and normal quark phases, respectively. This 
condition is satisfied at each point along the mid- 
dle part of the dark solid line (i.e., between □- and 
A-points). With increasing density (up to 2.56po)> 
the volume fraction of the hadronic phase decreases 
(down to 0.59). However, one does not reach a point 
where the fraction of hadronic phase would vanish 
completely. Instead, at baryon density about 2.56/?o 
(A-point in Fig. the mixed phase is replaced by 
another mixed phase which is made of the normal 
and 2SC quark components. At this point, the pos- 
itively charged hadronic component will be suddenly 
converted into a positively charged 2SC quark compo- 
nent. As a result of this rearrangement, the values of 
the baryon density and the energy density experience 
small jumps: the baryon density changes from about 
2.56po to 2.75po and the energy density changes from 
378 MeV/fm 3 to 415 MeV/fm 3 . After the transition, 
the fractions of the 2SC and normal quark phases are 
0.53 and 0.47, respectively. 

At higher densities, the mixture of the normal and 
2SC quark phases is the most favorable globally neu- 
tral construction in the model at hand. The volume 
fractions of the components in this mixed phase stay 
nearly constant with increasing density. 

The complete equation of state of hybrid baryon mat- 
ter is shown in Fig. [21 by solid line. For comparison, 
the equation of state for globally neutral quark mat- 
ter is also shown. As before, the points that indicate 
the beginning of two different mixed phases of hybrid 
matter are denoted by □ and A symbols. 

4 Star properties 

Here we use the equation of state of hybrid matter that 
is shown in Fig. [21 to construct non-rotating compact 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



e[GeV/fm 3 ] 

Figure 2. The equation of state for globally neutral hy- 
brid matter (solid line) and globally neutral quark matter 
(dashed line) . The points of the beginning of the two mixed 
phases are denoted by □ and A. 

stars. This is done by solving the well known Tolman- 
Oppenheimer-Volkoff (TOV) equations. 

The energy density profiles for hybrid stars with dif- 
ferent central energy densities are displayed in Fig. 01 
The star with the lowest central energy density in 
Fig- El e c = 210 MeV/fm 3 , is composed of pure con- 
fined hadronic matter, mainly neutrons, surrounded 
by a thin compact star crust consisting of leptons and 
nuclei. To get the equation of state for the crust, we 
use the results of Ref. [[201 for ps < 0.001 fm~ 3 and 
the results of Ref. [E] for 0.001 < p B < 0.08 fm~ 3 . 

The next energy density profile in Fig. corresponds 
to the central energy density e c = 370 MeV/fm 3 . We 
see that the corresponding star already has a rather 
large core (the radius is about 8 km) consisting of a 
mixture of hadronic and normal quark matter. The 
core of the star is surrounded by a layer of hadronic 
matter and a crust. 

In the model at hand, there are no stars with central 
energy densities in the window between 378 MeV/fm 3 
and 415 MeV/fm 3 . At e c > 415 MeV/fm 3 , a quark 
core (made of the mixed phase of the normal and 2SC 
components) forms at the center of the star. Two ex- 
amples of the corresponding energy density profiles are 
also shown in Fig. The star with the central energy 
density e c = 500 MeV/fm 3 contains a quark phase 
core with radius about 6 km. This core is separated 
from the layer of the hadron-quark mixed phase by 
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radius [km] 

Figure 3. Energy density profiles for hybrid stars. The 
locations of the interface between the two types of mixed 
phases are denoted by A, while the locations of the bound- 
ary between the pure hadronic phase and the hadron-quark 
mixed phase are denoted by □. 

a sharp interface (the corresponding point is denoted 
by A in Fig. EJ). The star with the largest possible 
mass within this model has the central energy density 
e c = 1392 MeV/fm 3 . This star is composed mainly of 
the quark core surrounded by a relatively thin layer 
of the hadron-quark mixed phase, as well as the pure 
hadronic phase and the crust on the outside. 

The results for the mass-radius relation of hybrid and 
quark stars are shown in Fig. 01 As one would ex- 
pect, the pure quark stars have much smaller radii and 
the value of their maximum mass is slightly smaller. 
The difference between hybrid and pure quark stars is 
mostly due to the low density part of the equation of 
state. This is also evident from the qualitative differ- 
ence in the dependence of the radius as a function of 
mass for the hybrid and quark stars with low masses. 
The corresponding hybrid stars are large because of 
a sizable low density hadronic layer, while the quark 
stars are small because they have no such layers. 

Our results for pure quark stars are comparable to 
those in Refs. [ E21 EH ES EH . Also, the maximum 
masses and the corresponding radii of the hybrid stars 
obtained here are similar to those of the strange hybrid 
stars of Ref. [ |^ , provided the strange quark mass 
is not very small (m s > 300 MeV) and the supercon- 
ducting gap is not too large (A < 50 MeV). At small 
values of the strange quark mass and/or large values 
of the superconducting gap, the strange hybrid stars 




2 1 < 1 < i < 1 < i < 1 < i < 1 < 1 

7 8 9 10 11 12 13 14 15 

R [km] 

Figure 4. The mass-radius relations for hybrid stars 
(solid line) and quark stars (dashed line). 

tend to have smaller maximum masses and smaller 
radii [EI]. 

5 Summary and outlook 

Here we constructed a realistic equation of state of 
non-strange baryonic matter that is globally neutral 
and satisfies the condition of /3-equilibrium. This 
equation of state might be valid up to densities of 
about 8po in the most optimistic scenario. In our 
construction, matter at low density (pb ^ 1.49po) is 
mostly made of neutrons with traces of protons and 
electrons. At intermediate densities (1.49po Pb ^ 
2.56po)j homogeneous hadronic matter is replaced by 
the mixed phase of positively charged hadronic mat- 
ter and deconfmed bubbles of negatively charged nor- 
mal quark matter. The volume fraction of normal 
quark matter gradually grows with increasing den- 
sity. Before the volume fractions of two phases become 
equal, the mixed phase undergoes a rearrangement in 
which the hadronic component of matter turns into 
the two- flavor color superconductor. At higher densi- 
ties (pb ^ 2.75po)j only the quark mixed phase exists. 
This latter is composed of about equal fractions of 
normal and 2SC quark matter. 

Previously, it was argued that 2SC quark matter could 
not appear in compact stars when the charge neutral- 
ity condition is imposed locally [ El • The main rea- 
son for this is the strong preference of the 2SC phase 
to remain positively charged. Our study shows, how- 
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ever, that the positively charged 2SC quark compo- 
nent appears naturally in a quark mixed phase at den- 
sities around 3po- The other component of the glob- 
ally neutral mixed phase is negatively charged normal 
quark matter. The corresponding construction turns 
out to be rather stable. In particular, we observe 
that the volume fractions of the two quark compo- 
nents remain approximately the same with changing 
the baryon density in a wide range. 

By making use of the equation of state of hybrid mat- 
ter, we construct non-rotating compact stars. We find 
that the largest mass hybrid star has the radius 10.86 
km, the mass 1.81M Q and the central baryon den- 
sity 7.58po- This star has a large (8 km) quark core, 
and a relatively thin outer layers of hadronic matter 
and a crust. One could speculate that the appearance 
of a large core is connected with the fact that the 
quark phase starts to develop at relatively low densi- 
ties, pb ~ 2.75poj in the model used. 

In the end, we would like to mention that perform- 
ing a systematic study of the model dependence of 
hybrid matter constructions with various color super- 
conducting phases would be an important task for fu- 
ture. Such a study will be crucial for resolving some 
apparent differences between the results currently ex- 
isting in the literature regarding non-strange as well 
as strange quark matter [ HI 1221 1251 121125] . 
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